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We study the low-energy properties of a triangular triple quantum dot connected to two non- 
interacting leads in a wide parameter range, using the numerical renormalization group (NRG). 
Various kinds of Kondo effects take place in this system depending on the electron filling iVtot, or 
the level position ed of the triple dot. The SU(4) Kondo behavior is seen in the half-filled case 
iVtot — 3.0 at the dip of the series conductance, and it causes a charge redistribution between the 
even and odd orbitals in the triangle. We show generally that the quasi-particle excitations from 
a local Fermi-liquid ground state acquire a channel symmetry at zero points of the two-terminal 
conductance, in the case the system has time-reversal and inversion symmetries. It causes the SU(4) 
behavior at low energies, while the orbital degeneracy in the triangle determines the high-energy 
behavior. At four-electron filling JV to t — 4.0, a local 5 = 1 moment emerges at high temperatures 
due to a Nagaoka ferromagnetic mechanism. It is fully screened by the electrons from the two 
conducting channels via a two-stage Kondo effect, which is caused by a difference in the charge 
distribution in the even and odd orbitals. 

PACS numbers: 72.10.-d, 72.10.Bg, 73.40.-c 



I. INTRODUCTION 

The Kondo effect in quantum dots is an active 
field of current research i 1 ' 2 ' 3 ' 4 It has been studied 
in various situations, such as Aharanov-Bohm (AB) 
rings,^' 7 -'^ tunneling junctions coupled to superconduct- 
ing hosts^ i 10 ' 11 ! 12 ! 13 ! 14 : 15 and a number of multi-dot and 
multi-orbital systemsi 16 ' 17 i 18 ' 19 Due to the variety of 
combinations of quantum dots that can be fabricated, 
and their high tunability, it is possible to study novel 
interplays between the strong electron correlation and 
quantum-mechanical interference effects. These were 
once studied as different topics in the field of the con- 
densed matter physics. 

The triangular triple quantum dot (TTQD) is an in- 
teresting system, which can demonstrate various types 
of the Kondo effects depending on an external field, 
the electron filling iVtot, and so on . 20 ' 21 For instance, 
a closed path along the triangle under a magnetic 
field can work as an AB interferometer^ Also in 
the case four electrons occupying the cluster of the 
triple dot, N tot = 4.0, the closed loop induces a lo- 
cal S = 1 moment due to a Nagaoka ferromagnetic 
mechanism^ Furthermore, the orbital degrees of free- 
dom cause a sharp conductance dip appearing at half- 
filling N tot = 3.0,2! which links to the SU(4) Kondo 
effects. These features distinguish the TTQD from a 
linear chain of three quantum dotS) 22 i 23 ' 24 i 25 ' 26 i 27 and 
from the other three-level systems ^21^^,32^3 E X p er j_ 
mentally, triple quantum-dot systems have been realized 
in AlGaAs/GaAs heterostructure ) 34 i 35 ' 36 self-assembled 
InAs^ 7 - and a single wall carbon nanotube^ The trian- 
gular trimmer of Cr atoms placed upon an Au surface is 
also a related systemi 39 ' 40 i 41 

The TTQD has a big parameter space to be explored, 



and the number of the conducting channels which are 
coupled to the triangle also gives interesting variations 
in low-energy properties. 42 ' 43 ^ 4 ' 45 In Refs. [43l and H the 
TTQD's connected, respectively, to two leads^ 3 - and one 
learli 4 - were studied. The main interest in these two re- 
cent works was the behavior at half-filling Ntot — 3.0, but 
the effects of the gate voltage ed which varies the electron 
filling in the TTQD were not examined. The other recent 
report Rcf. |45J addresses the effects of ed in the cases of 
the two and three leads, but the parameter region exam- 
ined was restricted in a small interaction U and a large 
dot-lead coupling T. In the previous works we studied the 
TTQD connected to two leads) 21 ' 42 and showed that the 
various types of the Kondo behavior clearly take place at 
different values of the gate voltage. The calculations were 
carried out, however, for a typical but just one parameter 
set which describes a large U and small T situation. 

The purpose of the present work is to study the fea- 
ture of the Kondo behavior in the TTQD connected to 
two leads in a wide parameter range. To this end, we 
calculate the conductance and scattering phase shifts at 
zero temperature, making use of the numerical renormal- 
ization group (NRG). Furthermore, we calculate also the 
TTQD contributions to the entropy and spin susceptibil- 
ity, and from their temperature dependence we deduce 
the Kondo energy scale. Our results reveal the precise 
features of the Kondo behavior in a wide parameter re- 
gion of the electron filling AT to t, interaction U, and the 
hybridization energy scale T. For instance, the 5 = 1 
moment is fully screened at low temperatures via two 
separate stages by the conduction electrons from the two 
non-interacting leads, which break the C^ v symmetry of 
the triangle. The two-stage screening process reflects the 
charge distribution in the even and odd orbitals which 
are classified according to the parity. It is confirmed 
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that the S = 1 Kondo behavior can be seen in a wide 
parameter region for N tot = 4.0. It is robust against the 
perturbations the typical energy of which is less than the 
finite energy separation between the Nagaoka and the 
first excited states for the triangular cluster. Specifically, 
the S = 1 Kondo behavior is not sensitive to the de- 
formations of the TTQD caused by the inhomogcncities 
in the level positions and the inter-dot hopping matrix 
elements. 

We show also that the fixed-point Hamiltonian, which 
describes the low-lying quasi-particle excitations from a 
local Fermi-liquid ground state, acquires an SU(2) sym- 
metry between the even and odd channels at zero points 
of the two-terminal conductance, in the case that the 
system has time-reversal and inversion symmetries. It 
means that if the phase shifts for the even and odd 
channels satisfy the condition S c — S Q — nn for n = 
0, ±1, ±2, . . ., then the low-energy behavior can be char- 
acterized by an SU(4) symmetric Fermi-liquid theory 
with the channel and spin symmetries, even if the orig- 
inal Hamiltonian does not have the SU(4) symmetry on 
a global energy scale. This happens quite generally for 
multi-dot and multi-orbital systems. In the case of the 
TTQD connected to two conducting channels, there are 
three zero points at most. Particularly, the one at half- 
filling iVtot = 3.0 shows a pronounced behavior as a sharp 
dip in the Kondo plateau. 21 ' 42 At this zero point the 
system has an SU(4) symmetry not only at low ener- 
gies but also at high energies, which is caused the or- 
bital degeneracy in the cluster of the regular triangu- 
lar. These symmetric properties cause the SU(4) Kondo 
behavio r 46 ' 47 ' 48 ' 49 seen in a wide temperature ranged 
Away from this zero point at the dip the NRG results 
of the thermodynamic quantities show that the channel 
symmetry is broken at low temperatures, and we observe 
the SU(2) Kondo behavior due to the spin. We find also 
that the sharp dip structure of the series conductance 
does not disappear from the Kondo plateau, even in the 
presence of infinitesimal deformations which break the 
Cs v symmetry, although the position of the dip shifts 
from the middle of the plateau as the deformation in- 
creases. 

We study also the characteristic temperature T*, at 
which a crossover to a singlet ground state take place, 
for a wide range of the electron filling iVtot- The result 
shows that T* becomes high at the conductance dip for 
Ntot = 3.0 due to the SU(4) Kondo behavior. It is en- 
hanced also at a mixed- valence point, where the ground- 
state energy for iVtot = 3.0 and that for iVtot = 4.0 
coincide. Furthermore, the screening temperature T* 
increases as the hybridization energy scale T increases, 
which would raise the experimental accessibility to the 
low-temperature Fermi-liquid region. 

The paper is organized as follows. In Sec. [UJ we gi ye 
the outline a description of the model and the formula- 
tion. In Sec. lIIIl we describe the low-energy channel sym- 
metry that the fixed-point Hamiltonian acquire at the 
zero point of the two-terminal conductance. The NRG 
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FIG. 1: Triangular triple quantum dot in (a) series and (b) 
parallel configurations. The dot labelled with i = 2 is referred 
to as the apex site, which has no direct connection to the leads. 



results for the ground-state properties are presented in 
Sec. IIV1 The results for temperature dependence of the 
TTQD contributions of entropy and spin susceptibility 
are presented in Sec. [V] The effects of the deformations 
which breaks the triangular symmetry are discussed in 
Sec. I VI I A summary is given in Sec. I VIII 



II. FORMULATION 

In this section we describe the model and the formu- 
lation which we have employed in the present work. We 
also describe some particular characteristics of the tri- 
angular triple dot in the two solvable limits, T = and 
U — 0, in order to see the essential ingredients that cause 
the various Kondo effects taking place in this system. 



A. Model 

We start with a three-site Hubbard model on a triangle 
that is connected to two non-interacting leads on the left 
(L) and right (R), at the sites labelled by i = 1 and 
i = No (=3), respectively, as illustrated in Fig. [TJ (a). 
The Hamiltonian is given by 



T~t — ^dot "F 1~(-dot F 'Hmix + ^lead j 

Kot=- E £%(4^+4<a 

<ij> a 

+££^4^* 



(f) 



i=l a 
No 



^dot = u e 



n d ,ta = d\ a d. 



(2) 



(3) 



n mix = J2{ v l A,„C La + v R J Noa C Ra + H.c. J , (4) 

Hlcad= E J2 e k C L* C kW ( 5 ) 



i/—L,R ka 
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Here, d\ a creates an electron with spin a at the i-th site 
in the dot, Cd,i the onsite energy in the dot, and U the 
Coulomb interaction. The hopping matrix elements ty 
between the dots are chosen to be positive (iy > 0). The 
conduction electrons near the dots, C va = ^2 k c kua /\f~N, 
can tunnel into the triangle via the hybridization term 
Eq. ([4|, which causes the level broadening r„ = ■npv 2 1 
with p the density of states of the leads. In the present 
work we consider the equal-coupling case T L = T R (= T), 
namely v L = v R (= v), and assume that the system has 
an inversion symmetry taking ed,i = td,3 (= £d), t\2 = ^23 
(= i) and ii3 (= t'). We shall refer to the dot at i = 2 in 
Fig. [1] as the apex site ed,2 = fapcx, and choose the Fermi 
energy Ep as the origin of the energy Ep = 0. 



B. Cluster of the triangular triple dot 

In order to see some characteristic features of the sys- 
tem, we hrst of all discuss a cluster of a triangular triple 
dot, described by H° ot + Ti^ ot - The eigenstates and en- 
ergies of this cluster give us a clue to understand the 
high-energy properties of the system. To be specific, we 
consider a regular triangle case, choosing e apcx = and 
t' = t. 

In the non-interacting case the one-particle states of 
H|J ot can be labelled by a wavenumber k (= 0, ±2n/3), 



and the eigenvalues are given by E 



(i) 



Note that AE^ > for U > 0, and in the two extreme 
limits at small and large U it takes the form 



AE^ 



U/3 
t , 



< [/< t 
0<t<^U 



(8) 



For weak repulsions, an infinitesimal U lifts the degener- 
acy of the singlet and triplet states in the non-interacting 
ground state. In the opposite limit, for large U the circu- 
lar motion along the triangle favors the magnetic 5 = 1 
ground state, and the energy separation AE^ is deter- 
mined by the hopping matrix clement t. If the cluster 
of the TTQD is connected to the leads, the local mag- 
netic moment will be screened by the conduction elec- 
trons showing the two-stage Kondo behavior described 
in Sec. EO 

As the level position ed decreases further, the occupa- 
tion number of the electrons increases. The ground-state 
energy Eq for a five-electron occupation and that for 
six electrons Eq are given by 



E, 



(5) 



E 



(4) 



E, 



(6) 



E 



(5) 



U + t + e d . 



(9) 



-2t COSfc + Erf, 



Specifically, at the point where the level in the TTQD 
takes the value of ed = —U — t, the ground-state energies 
for the three different fillings, TVtot = 4, 5 and 6, coincide. 
Thus the occupation number jumps from 7V to t = 4 to 6, 
as 6d crosses this value. 
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(i) 

k=0 



-2t + e d , 



E 



(i) 



= t 



(G) 



Phase shifts, conductances &: Friedel sum rule 



The excited states, for t > 0, are degenerate with re- 
spect to k = ±2tt/3. This degeneracy brings interesting 
varieties to the Kondo effects depending on the occupa- 
tion number iV tot , or e^, which can be controlled by the 
gate voltage. In the case that the number of electrons 
is iVtot = 3, namely at half- filling, the ground state for 
the triangular cluster has the four-fold degeneracy due to 
the orbital and spin degrees of freedoms. This is because 
the first two electrons occupy the one-particle state for 
k = 0, and then for the third electron there are four possi- 
ble ways to occupy one of the orbitals of k = ±27r/3 with 
the spin a =t, j- The four-fold degeneracy emerges also 
for a finite Coulomb interaction H^ ot , and could cause 
the SU(4) Kondo behavior. 

In the case that one additional electron is intro- 
duced into the triangular cluster, the ground state of 
H^ ot + H^ ot for a four-electron occupation iVtot = 4 be- 
comes a triplet with the total spin S = 1 for U > 
(see also appendix [5} . This is caused by the Nagaoka 
ferromagnetic mechanism^ and the energy separation 



A£W = E 



gj_ — -EgZj between the ground state E^L-^ 

.(4) 



U — 2t + 4td and the lowest singlet excited state EgJ_ is 
given by 



The charge transfer between the dots and leads makes 
the low-energy states of the whole system a local Fermi 
liquid, which can be described by renormalized quasi- 
particlcs. Particularly, in the case where the system has 
an inversion symmetry T L = T R , the ground-state prop- 
erties are characterized by the two phase shifts, 6 e and 
<5 , for the quasi-particles with the even and odd parities. 
At zero temperature, the series conductance g s in the 
two-channel configuration, which is shown in Fig. [1] (a), 
and the average number of electrons in the triple dots 
can be expressed, respectively, in the form of a Landaucr 
formula and the Friedel sum rulc j 23 ' 51 



2e 2 

g s = —r- sin 2 (S c - S ) , 
h 



(10) 



N **= E£<4a,> = -(^+So) . (ii) 



Furthermore, the parallel conductance g p for the current 
flowing along the horizontal direction in the four-terminal 
geometry, which is shown in Fig. [1] (b), can also be ob- 
tained from these two phase shifts 



AE^ = - 
2 



U + 3t- \/9t 2 + U 2 + 2Ut 



(7) 



2e 2 

g p = — (sin 2 S e + sin 2 5 Q ) . (12) 
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FIG. 2: (Color online) Conductances, phase shifts, and the 
filling JVtot of the triple dot in the non-interacting case U = 0.0 
are plotted vs ea/t for a regular triangle t = t, e ap cx = ea, 
and r/t = 0.12. 
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FIG. 3: Even and odd orbitals: the onsite potential of each 
orbital is given by e a o = eapex, e a i = ea—t', and en = td + t'- 



see in the upper panel (a) that the parallel conductance 
g p shows the peaks as the one-particle states defined in 
Eq. © cross the Fermi level, namely at the values of 



ed which satisfy E 



(i) 



0. We can see that the scries 



conductance g s has a zero point at ed = —t. This is 
due to the destructive interference, which reflects the 
value of the phase shifts S c = 3ir/2 and 6 = 7r/2 at 
this point. The overall ed dependence of the phase shifts 
2i5 e /7r, 2S /ir and local charge N tot defined in Eq. (jTT|) 
are shown in the lower panel (b). The phase shifts show 
a step of the height it as one resonance state crosses the 
Fermi level. Correspondingly the local charge also shows 
a staircase behavior. Note that two electrons with spin 
up and down enter the resonance state simultaneously in 
the non-interacting case. Among the three one-particle 
states k = and ±2ir/3, the lowest one with k = is an 
even-parity state, so that the even phase shift S e takes 
the first step at e<i ~ 2t. Then at td — —t electrons enter 
the other two one-particle states which can be classified 
into an even and odd states. 



D. Even and odd orbitals 



The phase shifts can be expressed in terms of the renor- 
malized matrix element iy t 22 i 23 i 24 (see also appendix [B]) 



^ + e«M*«+ReEi(0) , 



(13) 



where Sj(0) is the self-energy due to the interaction 
7^ ot . Specifically in the inversion symmetric case, the 
renormalized matrix element Uj for i,j (= 1,2,3) inside 
the TTQD takes the form 



{-Uj} 



e<2 



-t e. 



apex 1 

t e d 



(14) 



and the two phase shifts are given, respectively, by 



cot(5 P = 



■U-H- 2t 2 /e 



apex 



cot S = 



r 



(15) 



As a simple example, the conductances, the phase 
shifts, and the occupation number iVtot m the non- 
interacting case E^t = are plotted in Fig. [2] as func- 



tions of ed, for e a 



e d , f = t and T/t = 0.12. We 



The non-interacting Hamiltonian Ho = TL j ot + H m ix + 
7ii ca d can be written in a diagonal form, using the even- 
odd basis 



':>,a 



'Iff 



V2 



b 



'3a 



V2 



(16) 



and a, 



Oct 



L 2a- 



Here, for the linear combinations the 
labels and 1 are assigned in the way that is shown 
in Fig. [3] The couplings to the two leads by T breaks 
the C% v symmetry of the regular triangle, and lift the 
degeneracy of the even and odd states. In order to see 
this, we consider the spectral functions defined by 

M") = - - Im ( <(a 0a ; ajja, + «a 1(T ; a f lCT )) w ) , (17) 
A (u) = - - Im {{b la ;b\ a )) u , (18) 

where ((6;P)) U = -i J °° dt e l ( w + l0+ )* ({d(t),V}) is the 
retarded Green's function. These spectral functions for 
non-interacting electrons, U = 0, are plotted in Fig. 2] as 
functions of a; for T/t = 0.12 and = 0.0. The odd com- 
ponent A (u>) has a single peak at u> ~ t, and its width 
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is wider than that of the peak of A e (u>) emerging at the 
same position. This is because on the resonance at u ~ t 
the spectral weight of A e (cu) is dominated by the compo- 
nent of the apex site do, which is situated away from the 
leads, as shown in Fig. [3] In contrast, A (lo) represents 
the spectral weight of the level bi, which is directly con- 
nected to one of the leads. This geometry of the triangle 
and the two leads in the even-odd basis causes the dif- 
ference in the peak width. If the level position is lowered 
from e<2 = 0.0 which is the value chosen for Fig. [3J the 
peaks move towards the left side in the figure. As the two 
resonance states approach the Fermi energy uj = 0, the 
third electron enters the odd-resonance state first rather 
than the even one, because of the difference in the peak 
width. Then, both the resonance states become singly 
occupied when the two peaks just reach w = 0. Corre- 
spondingly, the phase shifts take the values 5 C = 3ir/2 
and 5 = tt/2, and the series conductance has a dip at 
e d /t = -l.Q in Pig. 

The interaction Hamiltonian defined in Eq. (J3j) can be 
rewritten, using the even-odd basis, in the form 



njV 
n dot 



Un a , ^ n a fli + — (n a ,it n a ,i[ + «mt n b,u) 
+ (^2 na l nb ' 1 ~ 2 S a ,i ■ Sb,l 
+ f ( at iT a uMiT + 6 1 A a U a lt) ' ( 19 ) 



Here, S a> i = E^' a L<^' a «T'A & tnc Pauli matrices, 
n a ,ia = a\ a a i(J1 and n a ,i = J2 a n a,ia- The operators n bM 
and Sb,i for the odd orbital b\ a are defined in the same 
way. We see in Eq. ([19)) that the interaction contains 
some long-range components between a\ and b\ orbitals. 
Specifically, the last two terms in the right hand side 
describe a charge transfer which breaks the charge con- 
servation in each of the even and odd channels. Thus for 
U 0, each of the phase shifts S e and S does not nec- 
essarily correspond, respectively, to (n a fl + n a- i) n/2 and 
( n b,i) tt/2. The Friedel sum rule given in Eq. (fTTj) holds 
in a way such that S c + 5 = (?i a ,o + n a,i + n b,i) 7r/2. 



III. LOW-ENERGY SU(4) SYMMETRY AT 
ZERO POINTS OF CONDUCTANCE 



In this section we show that the low-lying excitations 
acquire a special symmetry between channels at the zero 
points of the series conductance. This is a general prop- 
erty of a Fermi-liquid fixed point for a quantum-impurity 
system connected to two conduction channels. To be 
specific, we consider the case that the system has time- 
reversal and inversion symmetries. We use the discretized 
Hamiltonian Hm of NR G 52 ' 53 (see also appendix [Cj in 
this section not for a numerical purpose, but for show- 
ing explicitly the symmetric properties that emerge in a 
finite-size spectrum. 



o 

o 
c 

=3 

S_ 

-t— » 

o 
CD 

Q. 

CO 



9 p 
8 - 
7 - 
6 - 
5 - 
4 - 
3 - 
2 - 
1 - 
tu 
-3 



TT 

even 

odd 

U/(2rtt)=0.0 
r/t = 0.12 



1 I 1 1 1 

£, = 0.0 




-2 



■ 1 /f ° 
CO/t 



FIG. 4: (Color online) Spectral functions A e (ui) and A {ui) 
for non-interacting electrons U = 0.0, for td = 0.0, t' = t, 
eapex = £d, and F/t = 0.12. The even and odd components 
are defined in Eqs. (fIS jH |Tg )) . 



The two phase shifts S c and S for interacting electrons 
can be deduced from the fixed-point eigenvalues of the 
discretized Hamiltonian Hm defined in Eq. (jCip . follow- 
ing the procedure described in Ref. HiJ Conversely, for 
some particular values of the phase shifts, one can deduce 
that the system captures a special symmetry asymptot- 
ically at low energies. Particularly at the zero points 
of the series conductance, the low-lying energy spectrum 
for the even (s-wave) and odd (p-wave) channels become 
identical near the Fermi energy. It means that the quasi- 
particles of the local Fermi liquid acquire the channel 
degeneracy as well as the spin degeneracy. 

In order to see this, we go back briefly to a few es- 
sential points of the NRG. At a Fermi-liquid fixed point 
the low-lying energy spectrum of the many-body Hamil- 
tonian Hn defined in (jCljl can be reproduced by a dis- 
cretized version of free quasi-particlc Hamiltonian H^P , 



H iN) s A ( N - 1 )/ 2 



H -X- WW 
-"mix ' "lead 



cr=T,iM=e,o / 



(20) 



Here, the renormalized TTQD part, 7idot, in the first 
line is defined in Eq. (|B6j) . The second line is a diagonal 
form, for which /i (= e, o) is the channel index, and the 
label / is a quantum number for the eigenstates in each 
channel. The ground state for can be constructed 

in a such way that all the states for e^p < are filled by 



the 7; M<T particles, and the excitation energy of the 
quasi-particles near the Fermi level can be determined 
from the fixed-point eigenvalues of if at. 52 ' 53 

The phase shifts can be expressed as = — cot 5^ in 
terms of the parameter which determines the value 
of a Green's function at the Fermi energy, as shown in 
Eq. (|B8[) . The value of this parameter k p can be deduced 
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from the quasi-particle energy spectrum ^ 



2 -4± lim DA^grf) 



(21) 



at oatV c ;u 

Here, is a coefficient which is defined in (|C4[) . D is 
the half width of the conduction band, and g N {s) is a 
local Green's function at the interface n = of the dis- 
cretized version of the non-interacting lead ffffl given in 
(|C3j) . The equation (f2"Tj) represents the relation between 
the phase shift and the low-energy quasi-particle spec- 
trum in each channel \x. The right hand side of Eq. ([2Tj) 
converges to the same value of for any I, as long as 
is an excitation energy near the Fermi energy of the 
quasi-particles r ) l . Specifically, it is not necessary that 

e lu should be the lowest excitation energy, although the 
lowest one is preferable for a numerical purpose to make 
the corrections due to a residual interaction between the 
quasi-particles small^ 

One can deduce conversely the properties of the low- 
energy spectrum from the values of re e and k , using Eq. 
([2Tj) . Particularly at zero points of the series conduc- 
tance, the phase shifts satisfy the condition S c — S Q = mr 
for n = 0, ±1, ±2, . . .. This is equivalent to n c = K Q , and 
thus from Eq. (f2"Tj) it is deduced that the low-lying energy 
spectrum for the even and odd channels become identical 
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e} „ for all I near the Fermi 



for large N, namely 
level of the quasi-particles. Thus the quasi-particles ac 
quire a rotational symmetry in the channel (flavor) space 
asymptotically at low energies, in addition to the rota- 
tion symmetry of the spin. Therefore, has a SU(4) 
symmetry at the zero points of the series conductance. 

In the special case k c = and n a = 0, which cor- 
respond to S c = (n + l/2)n and S = {n 1 + 1 /2)7r for 
integer n and n', the low- lying energy states have a 
particle-hole symmetry in addition to the channel sym- 
metry, because the lead Green's function in the right 
hand side of Eq. (f2"Tj) is an odd function of the frequency 
g N (—e) = —g N (e), due to the particle-hole symmetry of 
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There are some other interesting cases in the unitary 
limit, at which the series conductance reaches g s = 2e 2 /h. 
If the two phase shifts take the values of 5 S = (n± 
and S = (n' =F l/4)ir, then k m becomes k c = =pl and 
k = ±1. Therefore in this case the low-energy excita- 
tions acquire a particle-hole symmetry between a parti- 
cle (hole) in the even channel and a hole (particle) in the 
odd channel, although there is no channel symmetry in 
this case. It can also be deduced from Eq. (f2"Tj) that for 
the other set of the phase shifts, 5 C = (n + l/2)w and 
8 = n'lr, particle-hole symmetry holds in each channel 
at low energies. 



IV. GROUND-STATE PROPERTIES 

We have carried out the NRG calculations in order 
to study the effects of the interaction on the transport 



and magnetic properties. Specifically we concentrate on 
the regular triangle case, for which ed = e ap cx and t' = 
t in the following. We will examine the effects of the 
deformations of the regular triangle caused by e ap ox and 
tf in Sec. ED 

The ratio of the inter-dot hopping matrix element t and 
the half width of the conduction band D are chosen typ- 
ically to be t/D = 0.1 in most of our calculations. The 
iterative diagonalization has been carried out by using 
the even-odd basis. For constructing the Hilbcrt space at 
each NRG step, instead of adding two orbitals from even 
and odd orbitals simultaneously, we add the one from the 
even part first and retain 3600 low-energy states after car- 
rying out the diagonalization of the Hamiltonian. Then, 
we add the other orbital from the odd part, and again 
keep the lowest 3600 eigenstates after the diagonaliza- 
tion. This truncation procedure preserves the inversion 
symmetry. We have chosen the discretization parameter 
to be A = 6.0. It has been confirmed that the conduc- 
tance for noninteracting electrons U = is reproduced 
sufficiently well with this calculation schemei 22 ' 23 i 24 



A. Results at a small F and large U 

The phase shifts and conductances for a relatively weak 
coupling T/t = 0.12 and strong interaction U/{2-Kt) = 1.0 
are plotted as functions of td/U in Fig.[5j The number of 
electrons iVtot m the triangle increases with decreasing ed- 
it shows plateaus for iVtot = 1, 2, 3 and 4, and then jumps 
to iVtot — 6.0 at ed — — 1.16E7 without taking a step for 
iVtot — 5.0. The same behavior has already been seen in 
the case of the triangular cluster, for which the ground- 
states energies for the three different fillings, iVtot = 4, 5 
and 6, coincide at ed = —t — U, as mentioned in Sec. Ill Bl 
The vertical dashed lines in (a) correspond to the values 
of ed at which the local charge varies discontinuously in 
the limit of V = 0. We can see that the staircase behavior 
of iVtot in this small T case is essentially determined by 
the one in this isolated limit. 

The triangular triple dot has the 5 = 1 local mo- 
ment due to the Nagaoka ferromagnetic mechanism at the 
plateau for AT to t — 4.0, as mentioned in Sec. Ill Al This 
moment is screened at low temperatures by the conduc- 
tion electrons from the two leads, and the ground state of 
the whole system becomes a singlet. We can see in Fig. 
[SJb) that the two conductances show a significant differ- 
ence at this filling N tot ~ 4.0. At -1.15J7 < e d < -0.85C/ 
the parallel conductance reaches the unitary limit of two 
conducting channels g p ~ 4e 2 //i, while the series conduc- 
tance g s is suppressed. We have seen a similar difference 
between g p and g s for U = near the peaks at ed = —t 
in Fig. H](a). We can see in Fig. O(a) that the two phase 
shifts are locked at the value S e ~ 37r/2 and 8 ~ tt/2 
in the wide region — 1.15J7 < ed < — 0.85J7, due to the 
Coulomb interaction. Therefore in this region, both the 
even and odd channels give a unitary-limit contribution 
to the four-terminal conductance g p , while the destruc- 
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FIG. 5: (Color online) Plots of NRG results as functions of 
e d /U for U/(2ixt) = 1.0, T/t = 0.12, t' = t and e apcx = e d . (a): 
the local charge JVtot, the phase shifts 28 e /ir and 2S /n. (b): 
the series (0) and parallel (■) conductances. The vertical 
dashed lines in (a) correspond to e^, at which JV to t in the 
limit of r = varies discontinuously. At the zero point of the 
series conductance for ~ — 0.59J7, the phase shifts take the 
values S ~ 7r/4 and <5 C ~ 5n/4. Inset in (b) shows the region 
near the zero point of g a . 



tive interference reduces the two terminal conductance 

9a- 

The difference between the two phase shifts 5 C — 5 Q is 
shown in Fig. [6l It crosses the value of ir at ed/U ~ 
— 1.15, —0.88 and —0.59. At these three points the series 
conductance becomes zero. The series conductance must 
have a maximum between two adjacent zero points, and 
thus a fine structure of <7 S seen at —1.2U < < — 1.1U 
can be regarded as one such example. It should be noted 
that in the non-interacting case only a single zero point 
appears as ed varies (see Fig. [2] (a)). 

In most of the other regions of the electron filling, the 
series and parallel conductances show a similar depen- 
dence. For instance, both the two conductances show a 
plateau with the value 2e 2 /h for 0.15J7 < < 0.3U in 



i — i — i — i — i — | — i — i — i — i- 

U/(2nt)=1.0 




0.5 0.0 
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FIG. 6: (Color online) Difference between the two phase shifts 
2(S C — 8 )/ir is plotted vs td/U for U/(2nt) = 1.0, for several 
values of T/t (□) 0.12 (0) 0.25, and (O) 0.5. At the points 
where the lines cross the horizontal line of 8 C — S c = 7T, the 
series conductance g a becomes zero. 



Fig. E^b). This Kondo behavior is caused by a S = 1/2 
moment of a single electron which occupies the lowest 



(i) 

fe=0' 



Then, for 



even-parity orbital with the energy E t 
— Q.3U < ed < 0.1U, two electrons fill the lowest or- 
bital, and both the series and parallel conductances are 
suppressed. As decreases further, for — 0.58U < < 
— 0.3f7, the third electron enters the odd orbital corre- 
sponding to the site &i which is illustrated in Fig.[3J The 
moment due to the third electron causes the SU(2) Kondo 
effect in this region. 

We can see in Fig. E{b) that there is a sharp dip in the 
series conductance at ed — —0.59U in the middle of the 
plateau for N to t — 3.0. This reflects a sudden change of 
the two phase shifts, seen in Fig. E{a). Near the series 
conductance dip, the difference between the two phase 
shifts <5 C — S Q varies from 7r/2 to 37r/2 continuously, taking 
the middle value of 7r at the zero point of g s . In contrast, 
the sum of the two phase shifts does not show a significant 
change, and is almost a constant 5 C + S Q ~ 37r/2 which 
corresponds to the occupation number A^ot — 3.0. Thus 
the parallel conductance varies very little in this narrow 
region, keeping the value g p ~ 2e 2 /h. Furthermore, the 
two phase shifts take the value 6 — tt/4 and 5 C ~ 57r/4 
at the zero point of g s in the center of the dip. These 
values of the phase shifts are consistent with the one re- 
ported for a capacitivcly coupled double dot with a single 
electron in an SU(4) symmetric limit 4^ As mentioned in 
Sec. IIII1 the quasi-particle excitations have channel de- 
generacy as well as the spin degeneracy generally at zero 
points of the series conductance, and the system has an 
SU(4) symmetry of a Fermi-liquid origin at low energies. 
In the case of the TTQD, the system also has an SU(4) 
symmetry of a different origin at high energies, which is 
caused by the four-fold degenerate eigenstates of the iso- 
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FIG. 7: (Color online) Filling of even and odd orbitals, (n a ,o), 
(n a .i) and (n^i) are plotted vs ea/U for U/(2nt) = 1.0 and 
r/f = 0.12. The dashed lines correspond to the phase shifts 
2S c /tt and 2S /n. Just at the zero point, ed ~ — 0.59J7, of 
the series conductance in the middle of the crossover region 
the phase shifts take the values 8 ~ t/4, <5 c ~ 57r/4, and 
three electrons distribute in these orbitals as (n a ,o) — 1-0, 
(n a ,i) ~ 1.28, and (n M ) ~ 0.72. 



lated triangular cluster with three electrons, mentioned 
in Sec. Ill Bt We will see in Sec. W\ from the results of the 
thermodynamic quantities that the four-fold degenerate 
local degrees of freedom remain free at high tempera- 
tures, and are screened at low temperatures by the con- 
duction electrons keeping the SU(4) symmetric behavior. 

The sudden change of the two phase shifts also implies 
that a charge redistribution takes place in the TTQD. 
One might expect from the behavior of S e and S that 
the third electron moves away from the odd orbital to 
enter the second even orbital for — Q.8U < < —0.6U. 
Each of the two phase shifts, however, does not necessar- 
ily correspond to the occupation number of the orbitals 
of each parity, N cvcn = (n a ,o + n a ,i) or N odd = (n^i ), in 
the interacting case U ^ 0, as mentioned in Sec. Ill Dl For 
this reason, we have calculated also the expectation value 
of (n a ,i) and (n^i) with the NRG in order to clarify pre- 
cisely the charge distribution. The results are presented 
in Fig. [7l We see that an intuitive estimate, which as- 
sumes iVeven ~ 2S c /tt and 7V dd ~ 2S /tt, approximately 
works for < — 0.85L/ and — 0.6U < e<j. For instance, 
the third electron really enters into the odd orbital to 
take the occupation number of the value, (nb t i) ~ 1.0, 
for -0.58C7 < e d < -0AU on the right side of the se- 
ries conductance dip. Such a separation of the Friedel 
sum rule, however, does not make sense in the central re- 
gion. We can see for —0.75 < td/U < —0.6 that the local 
charges are redistributed between the even and odd or- 
bitals, and take the values (n a ,o) — 1-0, (n a ,i) — 1-5, and 



phase shifts take the values 25 c /ir ~ 3.0 and 28 /ir ~ 0.0. 
Furthermore, at the zero point of the series conductance 
the phase shifts take the values S ~ ir/4 and 6 Q ~ 57r/4 
as mentioned, while the occupation number for the even 
and odd orbitals become (re o) — 1-0, (n a ,i) — 1-28, and 
(n M ) ~ 0.72. 

The charge distribution in real space can also be de- 
duced, using the operator identities n<j,2 = n a,o and 
nd,i + n d.3 = n a ,i + ^6,i which follow from the defini- 
tion in lit 



(nd,i) = (nd.3) = (n a ,i +n&,i}/2 



(22) 



0.5. Thus N c 



2.5 and A^dd — 0.5, while the 



We have confirmed that the numerical value of (n a ,o) 
and (n a .i + n.fe 1 i)/2 arc very close in the whole region 
of e d plotted in Fig. [7] Thus the occupation number of 
the apex site, (71^2), is not much different from (nd,i), 
or (fid 3). Consequently, for small T, the charge distri- 
bution is almost homogeneous in real space. The charge 
inhomogeneity in the even-odd basis, however, affects the 
transport and magnetic properties very much. This is be- 
cause the orbitals a± and b\ are connected directly to the 
leads as shown in Fig. [31 while the apex site does not 
have a direct connection to the leads. 

We can see also at -1.15Z7 < e d < -0.83J7 in Fig. 
that four electrons in the TTQD are distributed in the 
even and odd orbitals such that (n a> o) — 1.3, (n a ,i) — 1.7, 
and (715,1) — 1-0. These values are close, respectively, to 
4/3, 5/3 and 1, which are the occupation numbers in the 
cluster limit V = (see appendix 



B. Large V cases 

The conductances and phase shifts discussed in the 
above are the results obtained for a relatively small hy- 
bridization T/t = 0.12. For this reason the structures 
of the plateaus and valleys can be seen very clearly. We 
next examine how the strength of the couplings T be- 
tween the triangle and two leads affects the ground-state 
properties. For comparison with the data shown in Fig. 
[SJ we have carried out the calculations for several larger 
values of r keeping the interaction U the same value of 
U/(2nt) = 1.0. The results arc shown in Fig. [5J The 
couplings are chosen such that T/t = 0.25 for (a) and 
(b), and T/t = 0.5 for (c) and (d). 

As a general trend, the structures seen in the phase 
shifts and conductances are smeared gradually as the 
couplings get stronger. We see in the upper panels (a) 
and (b) of Fig. [51 however, that the structures are still 
distinguishable by eye in the case of T/t = 0.25, and 
the feature of each line can be compared with the corre- 
sponding line in Fig. [SJ Specifically, for A^ot — 3.0 the 
broadening makes the dip of the series conductance near 
€d — — 0-6U wider. This structure reflects the behavior of 
the phase shifts at -0.8U < e d < -0.6U in Fig.[8j(a). In 
this region the odd phase shift S Q still shows a clear con- 
caving valley, although it becomes shallower than that 
in the case of Fig. [5] (a). This valley is smeared further 
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FIG. 8: (Color online) Phase shifts and conductances for large 
hybridizations (upper panels) T/t = 0.25 and (lower panels) 
T/t = 0.5; U/(2irt) = 1.0 is fixed. The vertical dashed lines 
in (a) and (c) correspond to ed, at which JVtot in the limit of 
T = varies discontinuously. 



FIG. 9: (Color online) Phase shifts and conductances for small 
values of U (upper panels) U j (2irt) = 0.5 and (lower panels) 
U/{2wt) = 0.2; T/t = 0.12 is fixed. The vertical dashed lines 
in (a) and (c) correspond to ed, at which JVtot for T = varies 
discontinuously. 



in the case of T/t = 0.5, and the fine structure of g s at 
-0.8C/ < e d < -0.6J7 is not visible in Fig. M (d). Par- 
ticularly, the conductance plateau on the left side of the 
dip almost vanishes. We can see, nevertheless, in Fig. [6] 
that the difference between the two phase shifts S c — S 
for T/t = 0.5 crosses the value of ir at ea/U ~ —0.73, 
—0.82, and —1.15. It means that there still exist three 
zero points of the series conductance, and g s has a weak 
maximum between the two adjacent zero points. 

In contrast to the series conductance <? s , the parallel 
conductance g p is less sensitive to T. Particularly, we can 
see that in Fig.[8]that the broad peaks of g p at iVtot — 4.0 
and JVtot — 1-0 keep their essential features seen in the 
small r case in Fig. [SJ This is because the energy gain 
of the Nagaoka state AE^ 4 \ which is defined in Eq. (J7J), 
is of the order of the inter-dot hopping matrix element 
t for large U. It takes a value of AE^ ft = 0.735 for 
U/(2nt) = 1.0, and thus in the present case AE^ 4 ' /t is 
still much larger than the energy scale of the hybridiza- 
tion T/t = 0.5. 

We can also see another quantitative change in Fig. [H] 
(d) that the difference between g p and g s at N tot ~ 1.0 
becomes visible for large T. A similar tendency was seen 
also for the conductances of a linear triple dot^ 



C. Small U cases 

We next consider the ground-state properties for weak 
interactions. For comparison with the case of Fig. [5l we 
have carried out the calculations, taking U smaller than 
U/(2irt) = 1.0 and choosing the same value of T/t = 0.12 
for the hybridization. The results are plotted in Fig. O 
the upper panels (a)-(b) for U /(2irt) = 0.5, and the lower 



panels (c)-(d) for U/(2-Kt) = 0.2. The vertical dashed 
lines in the panels (a) and (c), show the values of the 
ed at which the occupation number iVtot changes in the 
isolated limit T = 0. 

We can see that the area for three-electron filling be- 
comes broad as U increases. Note that in the non- 
interacting case the occupation number AT to t does not 
take the steps at N to t — 1-0, 3.0 and 4.0 as shown in 
Fig. H£b). Therefore the effects of the Coulomb interac- 
tion are particularly important for these electron fillings. 
The conductance peak at JVtot — 1-0 becomes wider as 
U increases, and the peak shape deviates from a simple 
Lorentzian form which is seen in Fig. [2](b) for U = 0. This 
is caused by the SU(2) Kondo effects due to the moment 
of a single electron occupying the local one-particle state 
^i=o m ^ ne triangle, as mentioned. 

We see in (a) and (b) of Fig. [9] that the overall feature 
of the phase shifts and the conductances for U / (2irt) = 
0.5 are similar to the those for U/{2irt) = 1.0 shown 
in Fig. \5\ For instance, the dip structure of the series 
conductance is clearly seen for JV to t — 3.0 in (b), which 
reflects the behavior of the phase shifts at —0.8 < ed/U < 
—0.6 in Fig. [9] (a). Particularly, the peak of g s on the 
left side of the dip is caused by a concaving behavior 
of 5 . The occupation number JV to t; however, no longer 
shows a clear staircase behavior at three-electron filling 
for U/ (2irt) = 0.5. Correspondingly, we see in Fig. [5] (b) 
that the parallel conductance g p has a shoulder behavior, 
instead of a well-defined plateau, at —0.8 < ed/U < —0.6. 

These plateau and valley structures almost disappear 
for the smaller interaction U/(2-Kt) = 0.2, as shown in 
the lower panels (c) and (d) of Fig. [51 Specifically, 
the region for three-electron filling, which corresponds 
to —1.2 < ed/U < —1-1, becomes very narrow, as men- 
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tioncd. Nevertheless, we have confirmed that the differ- 
ence between the two phase shifts 5 C — S crosses the value 
of 7r at the points ed/U ~ —1-2, —1.4 and —1.7, which 
correspond to the zero points of g s . The series conduc- 
tance has a weak peak between the two adjacent zero 
points, although it is not visible on the scale of Fig. [9j 

The parallel conductance at four-electron filling, which 
is shown in (b) and (d) of Fig. [5J keeps the features of 
the 5 = 1 Kondo behavior seen for larger U in Fig. [5] 
This is because the Coulomb interaction is still larger 
than the inter-dot hopping matrix element t, even for 
U/(2irt) = 0.2. Correspondingly, we can see in Fig. [5] (c) 
that the even phase shift shows a clear plateau behavior 
for S e ~ 37r/2 at -1.7 < e d /U < -1.3, although the odd 
phase shift S a behaves moderately. Note that the energy 
gain of the Nagaoka state is given by AE^ = 0.319 1 for 
C//(27rf) = 0.2, and AE^ = 0.563 1 for U/(2nt) = 0.5. 
Therefore, in both of these two cases the energy AE^ is 
greater than T (= 0.12t), and it makes the 5=1 Kondo 
behavior robust. 



V. THERMODYNAMIC QUANTITIES 

We present the NRG results for the magnetic suscepti- 
bility and the entropy of the TTQD in this section. The 
Kondo behavior in a wide energy scale can be investi- 
gated from the temperature dependence of these ther- 
modynamic quantities. Particularly, we study how the 
screening of the local moment is completed in the case 
of the 5=1 high-spin Nagaoka state at AT to t — 4.0. We 
investigate also the SU(4) Kondo behavior, which takes 
place at the series conductance dip in the middle of the 
plateau at AT to t — 3.0. 



A. Entropy &z spin susceptibility 

The free energy of the whole system consisting of the 
dots and leads is given by F = — Tlog [Tre - W / T l. The 
contribution of the quantum dots on this free energy can 
be extracted, by subtracting the contribution of the non- 
interacting leads, Ficad = —Tlog [Tr e _Wlcad / T ] , from the 
total free energy^ 
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FIG. 10: (Color online) Entropy S and spin susceptibility 4Tx 
in the middle of the charge step for iVtot — 4.0 are plotted vs 
bg(T/D) for t' = t, eapox = e d and F/t = 0.12, where D 
the half-width of the conduction band. In (a) both S and 
4Tx are plotted for U/(2-irt) = 1.0 and e d = -1.01(7. In 
(b) entropy is plotted for several values of U : (dash-dot line) 
U/(2iYt) = 0.2, e d = -1.52(7, (dotted line) U/(2ivt) = 0.5, 
e d = -1.01(7, (dashed line) U/(2irt) = 1.0, e d = -1.01(7, and 
(solid line) U/(2ivt) = 1.5, e d = -0.99(7. 



Hamiltonian H defined in Eq. ([T]), 



AF = F - Fi 



lead 



(23) 



Then, the entropy 5 and the spin susceptibility x due to 
the dots can be obtained from 



S = - 



dAF 

dT ' 



\ 



d 2 AF 
dH 2 



(24) 



Here, the magnetic field H is introduced both for the 
dots and leads, adding the Zeeman energy 7i ex to the 



"Hex = - {nd,i\ - n d ,ii) H 

i=l 

u=L,R k 



These thermodynamic quantities can be evaluated us- 
ing the eigenvalues of the discretized Hamiltonian of the 
NRG £2 
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B. Screening of S — 1 Nagaoka high spin 

The entropy S and 4%T are plotted in Fig. [TO] (a) as 
functions of log(T/D) for a parameter set. U / (2irt) = 1.0 
and e d = — l.Olf/, which corresponds to the middle of 
the charge step for N to t — 4.0 in Fig. O We can see 
that the screening of the spin S = 1 local moment is 
achieved via two separate stages. The first and second 
stages take place at T/D ~ 10" 8 and T/D ~ 10~ 29 , re- 
spectively. At high temperatures 10~ 7 < T/D < 10~ 3 , 
the entropy shows a plateau of the value of S ~ log 3, 
which represents the contribution of a free unscreened 
S = 1 moment. Then, at intermediate temperatures 
10 -27 < jy£j < iq-10 ) tne entropy takes the other value 
S ~ log 2 corresponding to a spin 1/2 moment. In this 
region, the spin susceptibility also shows a Curie behavior 
4xT ~ 1.0 =S> 4s(s + l)/3 with a coefficient s = 1/2 due 
to a free spin. Therefore in the intermediate temperature 
region, half of the 5 = 1 moment is screened by the con- 
duction electrons from one of the channel degrees of free- 
dom, and the local moment is still in an under-screened 
situation with the remaining free spin of s = 1/2. The full 
screening is completed at low temperature of the order 
T/D ~ 10~ 30 in the case of Fig. [TO] (a). This crossover 
energy scale to the low-energy Fermi-liquid regime can 
be regarded as a Kondo temperature for the 5 = 1 local 
moment in the present case. We refer to the crossover 
energy scale to the singlet ground state as the screening 
temperature T* , and use this definition for a wide range 
of electron filling < N to t < 6. 

The reason that the two stages are required to reach 
the Kondo singlet state can be understood from the dis- 
tribution of the 5 = 1 moment in the triangular triple 
dot. Specifically in the limit of T = 0, the partial moment 
m a /b.ii which constitutes the 5=1 moment as shown in 
appendix El distributes in the even-odd orbitals ao, a\ 
and bi with the weights 1/3, 1/6, and 1/2, respectively. 
Among the three orbitals, a\ and b\ are coupled directly 
to the lead, as illustrated in Fig. [3] Therefore, the partial 
moments in these two orbitals can be screened directly 
by the conduction electrons at the first stage. In con- 
trast, the other orbital ao has no direct tunneling matrix 
element to the leads, and thus the partial moment in the 
ao orbital has to be screened indirectly by the conduc- 
tion electrons coming over the a\ orbital which is already 
filled up to (n a ,i) — 1.7, as mentioned in IIV Al This 
amount of the occupation number is caused mainly by 
the difference between the onsite potential e a o = and 
e a \ = 64 — t in the even-odd basis, namely e a \ is deeper 
than e a o (see also appendix IA"]) . Therefore, in this sit- 
uation the charge and spin fluctuations must be nearly 
frozen in the a\ orbital, and the screening temperature 
for the partial moment in the apex site ao becomes low. 
These features of the screening at the second stage are 
analogous to those of a super-exchange mechanism. 

We see in Fig. [TO] (a) that the essential feature of the 
screening can be observed either in S or T\- At the 
transient regions, however, S varies more clearly than 



T\ does. Therefore, in order to see how the Coulomb in- 
teraction U affects the two-stage screening of the 5 = 1 
moment at N to t — 4.0, the entropy S is plotted in Fig. 
[TO1 (b) for several values of U, choosing the parameter 
set (U/(2nt), ea/U) in a way such that it corresponds 
to the middle of the plateau of the parallel conductance, 
solid line: (1.5, -0.99), dashed line: (1.0, -1.01), dot- 
ted line: (0.5, -1.01), and dash-dot line: (0.2, -1.52). 
The hybridization is chosen to be Y/t = 0.12. Note that 
the ground-state properties for these two small U cases 
were presented in Fig. [5] We see in Fig. [TO] (b) that the 
crossover temperature decreases as U increases. Partic- 
ularly, the second screening stage is more sensitive to U 
than the first stage which occurs at higher temperature. 
We can also see that the crossover temperature of the first 
stage and that of the second one become close for small 
U. For instance, the two stages take place almost suc- 
cessively for U/(2wt) = 0.2. Note that the value of the 
coupling between the TTQD and leads chosen for Fig. 
[TO] is relatively small Y/t = 0.12, so that the screening 
temperature to reach the singlet state becomes very low 
in the present case, especially for large U. The screen- 
ing temperature, however, rises as Y increases. It would 
make T* an accessible value in experiments. We will ex- 
amine the r dependence of T* in Sec. IV Dl 



C. SU(4) Kondo effect at the dip for 7V to t = 3.0 

In this subsection we present the results of the ther- 
modynamic quantities at three-electron filling A to t — 3.0. 
In Fig. 1 1 1 1 the entropy and spin susceptibility are plot- 
ted vs log(T/£>), for U/{2nt) = 1.0, and Y/t = 0.12. 
The impurity level is chosen in a way such that (a) 
ed = — 0.5897J7 at the zero point of g s in the middle 
of the sharp dip seen in Fig. [5] and (b) a = — 0.65C/ 
(— 0.55J7) at a point on the left (right) of the dip. 

We can see in Fig. [TTJ (a) that at high temperatures 
10~ 5 < T/D < 10~ 2 the entropy is a constant S ~ log 4, 
which is caused by the four-fold degeneracy of the lo- 
cal states in the isolated TTQD cluster mentioned in 
Sec. IIIBI Furthermore, the spin susceptibility shows the 
Curie behavior \ — s ( s + l)/(3r) with s = 1/2 at high 
temperatures. Then, at low temperatures both S and \T 
decrease with T. Specifically, the screening of the local 
four-fold degenerate states is completed in a single stage 
at T* /D ~ 10~'. This behavior is observed just at the 
zero point, at which the two phase shifts arc 5 C ~ 57r/4 
and S a ~ 7r/4. At the zero point the low-energy excita- 
tions from the singlet ground state asymptotically have 
the channel symmetry in addition to the global SU(2) 
symmetry of the spin, as shown in Sec. lIIII Therefore, the 
low-temperature behavior of these thermodynamic quan- 
tities must reflect the low-energy SU(4) symmetry of the 
fixed-point Hamiltonian. We refer to this zero point of 
the series conductance at N to t ~ 3.0 as the SU(4) sym- 
metric point in the following. It should be noted that T* 
is enhanced at this point due to the channel degeneracy 
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FIG. 11: (Color online) Entropy S and spin susceptibility 
4Tx in the step for JV tot ~ 3.0 are plotted vs log(T/D) for 
t' — t, e ap cx = £d and T/t — 0.12. Impurity level is chosen to 
be (a) just at the series conductance dip = —0.5897(7, (b) 
on the left = —0.6517 and the right = —0.55(7 sides of 
the dip. 



For the values of td away from the dip, the entropy 
shows the two-stage behavior. We can see in Fig. [TT] (b) 
that the entropy still shows the value of S ~ log 4 at 
high temperatures even in the case away from the dip. 
At intermediate temperatures, however, it takes another 
plateau with S ~ log 2, and the spin susceptibility shows 
the Curie behavior with the coefficient of s = 1/2. There- 
fore, the orbital part of the degrees of freedom is frozen at 
intermediate temperatures, and the spin degrees of free- 
dom still remains free. Furthermore, away from dip the 
low-lying energy excited states no longer have the chan- 
nel symmetry. The full screening of the remaining spin 
1/2 is completed at T* / D ~ 10~ 21 for e d = -0.65(7, and 
T*/D ~ 10~ 8 for ed = —0.55(7. These temperatures are 
lower than the screening temperature at the SU(4) sym- 
metric point, but still much higher than T* for the S = 1 
Kondo effect at N to t = 4, for the same parameter values 
of U, t, and T. 

Furthermore, the screening temperature on the right- 
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FIG. 12: (Color online) Entropy vs log(T/D) for several val- 
ues of V; U/(2irt) = 1.0 is fixed. Impurity level is chosen 
to be (a) Ed/U = —1.01 which corresponds to iVtot — 4.0, 
and (b) at the dip in the middle of the step for JVtot — 3.0 
(e d /U = -0.55, -0.62, and -0.73 for T/t = 0.12, 0.25, and 
0.5, respectively). 



hand side of the SU(4) symmetric point is higher than 
that on the left-hand side. The difference reflects the 
charge redistribution, taking place at the conductance 
dip as discussed in Sec. IIV Al As we saw in Fig. [71 
the distribution is almost homogenous on the right side; 
(n a fi) — ( n a,i) — (nb.i) — 1-0. In contrast, it becomes 
inhomogeneous on the left side; (n a , ) ~ 1.0, (n a ,i) — 1.5 
and (nfc.i) — 0.5. In the second screening stage the par- 
tial moment in the ao orbital, which is singly occupied 
on average, has to be screened by the electrons coming 
over the ai orbital. Therefore, T* must be affected by 
the filling of the ai orbital. Specifically, the charge and 
spin fluctuations at the a\ orbital, which is necessary for 
the second stage screening, tends to be frozen as (n 0) i) 
approaches to 2.0. The situation is similar to that in 
the S = 1 Kondo case, mentioned in Sec. IV Bl We will 
discuss this point further in the next subsection. 

The conductances which we discussed in Sec. IIVI arc 
the results obtained at zero temperature. If the tem- 
perature is raised, the series conductance at the SU(4) 
symmetric point will increase. Simultaneously, the series 
conductance plateaus on the left- and right-hand sides 
will decrease. Then, the dip structure will disappear at 
the temperature of T* on the left side, because it is the 
lowest energy scale near the symmetric point. Note that 
the screening temperature T* itself varies depending on 
the parameters T, t and U. In order to preserve the typ- 
ical structure of the dip, however, there are upper and 
lower bounds for T/t and U/t, respectively, as we can see 
in Figs. |U and H 



D. Dependence of the screening temperature on T 

We consider in this subsection the dependence of the 
screening temperature T* on the hybridization T and the 
level position e^. For comparisons, we have calculated 
the entropy for several values of the hybridization T/t = 
0.12, 0.25 and 0.5, for a fixed value of the interaction 
U/(2irt) = 1.0. The ground state properties for these 
parameter sets were presented in Figs. [S] and [5] 

Figure [T2l shows the results obtained at (a) the paral- 
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FIG. 13: (Color online) The crossover energy scale T* as a 
function of ea for several values of T/t = 0.12,0.25, and 0.5; 
U I (2-Kt) = 1.0 is fixed. T* is defined as the temperature in the 
middle of the low-energy crossover observed in the entropy. 
The vertical dashed lines correspond to e^, at which TVtot in 
the limit of F = varies discontinuously. 



lei conductance plateau for AT tot ~ 4.0 (ed/U = —1.01), 
and (b) the zero point of the series conductance for 
iVtot — 3.0. Therefore, we can see in (a) the T dependence 
of the two-stage screening of the 5 = 1 moment of the 
Nagaoka state, and in (b) the single-stage SU(4) Kondo 
behavior. In both of these two cases, the entropy for the 
intermediate coupling T/t = 0.25 retains the structures 
which we observed in the weak coupling case T/t = 0.12. 
For large hybridization T/t = 0.5, however, the charge 
transfer between the leads and the TTQD brings the sys- 
tem into a mixed- valence regime, and the typical struc- 
tures are smeared. Furthermore, we can see in (a) and 
(b) that the crossover temperature T* is sensitive to the 
value of the hybridization, and it increases rapidly with 
the coupling T. 

We have calculated the entropy also for different fill- 
ings, and determined the screening temperature T* from 
the mid-point of a low-energy crossover to a singlet 
ground state. The energy scale T* defined in this way 
can be regarded as a Kondo temperature in the case that 
the system has a well defined local moment at high en- 
ergies, which is the case of TV-tot — 1.0, 3.0, and 4.0. The 
results of T* are plotted as a function of td in Fig. [T3J We 
see that at the single-electron filling for 0. 1 < e d < 0.3/7 
the screening temperature is larger than that of the other 
regions of the filling. This is because at this filling the 
single electron enters mainly the a\ orbital which is con- 
nected directly to one of the leads, as illustrated in Fig. 

13 

In the region of the three-electron filling N to t — 3.0, the 
screening temperature is enhanced at the SU(4) symmet- 
ric point, td — —0.59/7, where the series conductance has 
a dip. We see in Fig.[I3]that T* in the weak hybridization 
case T/t — 0.12 has a local minimum at e d ~ —0.65/7 on 
the left side of the symmetric point, and T* becomes very 



low near the minimum. This is because the fraction of the 
local spin in the ao orbital must be screened indirectly by 
the conduction electrons coming through the ci\ orbital 
which is filled already up to (n a ,i) — 1-5. As we see in 
Fig. [3 there is a significant difference between the charge 
distribution on the left and right sides of the symmetric 
point. On the right side, at -0.58 < e d /U < -0.4, the 
local moment is mainly in the odd-parity b\ orbital which 
has a direct tunneling matrix element to one of the leads, 
as discussed in Sec. IIV Al For this reason, the screening 
temperature on the right side of the dip becomes larger 
than that on the left side. 

We see also in Fig. [T2]that the screening temperature 
is enhanced significantly at td — —0.82/7 for T/t = 0.12. 
This seems to be caused by the charge fluctuations be- 
tween the two different fillings of N to t = 3 and AT tot = 4. 
We note that near this point there is another zero point of 
the series conductance at e d /U ~ —0.88, as we see in Fig. 
[5] At four-electron filling, T* becomes small again par- 
ticularly at ed — —0.95/7. This is because in the second 
screening stage of the S = 1 moment the partial moment 
in the ao site needs to be screened by the electrons from 
the leads which have no direct hopping matrix elements 
to the apex site, as mentioned in Sec. IV Bl We see also 
that the screening temperature depends sensitively on T, 
as well as on U. Therefore, in order to observe the two- 
stage behavior experimentally, T should be tuned to a 
large value, but still has to be T < max(Z7/3, t) in order 
to retain the typical structures. This criterion is deter- 
mined by the energy scale AE^ ~ max(/7/3,t), defined 
in Eq. ©. 



VI. DEFORMATIONS OF TRIANGLE 

So far, we have assumed the full C% v symmetry for the 
TTQD. Although it has already been broken partly by 
the connection to the two leads, the perturbations which 
break the triangular symmetry should be examined fur- 
ther, since real TTQD systems must be deformed from 
the regular structure to a certain extent. We consider in 
this section how an asymmetry caused by the inhomo- 
geneity in the level positions and the inter-dot hopping 
matrix elements affect the ground-state properties. 

In the presence of a deformation which is described 
by e a pex and t', the one-particle energies of the isolated 
TTQD vary from those for the regular triangle given in 
Eq. (JHJ) to the following forms 



t-t(I) _ e apex + £d — t , 
^c,± - n ± 



"apex 



- e d + 1> 



E^ = e d + t' 



2t 2 , 

(26) 
(27) 



Here, E^l and E ( 1] are the eigenvalues for the states 
with even and odd parities, respectively. Note that the 
system still has the inversion symmetry in the present 
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FIG. 14: (Color online) Phase shifts and conductances are 
plotted vs ed, keeping the difference between e apcx and 
unchanged; (upper panels) e apcx — £d = —0.05(7 and (lower 
panels) e apox - e d = -0.1(7, with U/(2irt) = 1.0, T/t = 0.12 
and t' = t. 




FIG. 15: (Color online) Phase shifts and conductances are 
plotted vs ed, keeping the difference between e apcx and ed 
unchanged; (upper panels) e apcx — £d = 0.05(7 and (lower 
panels) e apcx - e d = 0.1(7, with (7/(2vrf) = 1.0, T/t = 0.12 
and t' = t. 



case. Among these three one-particle states, E c _ be- 
comes the lowest energy for small perturbations. The 
other two excited states become degenerate, E^ 1 } = Eo , 
in the regular triangle case with t' = t and e ap ex = ^d- 
The deformations lift this degeneracy, and the first order 
correction with respect to Se = e apox — ed and 5t = t' — t 
is given by 

*\**-\*t- (28) 

Thus, the energy of the even state E^i. becomes higher 

(lower) than that of the odd one Eo for Se — 2St > 
(< 0). It reflects the fact that the apex site with e ap cx 
belongs to the even part of the basis. Furthermore, t' 
shifts the onsite potential of the b\ orbital upwards, and 
that of the a\ downwards. In the following, we examine 
effects of Se and St separately, in a wide range of the 
electron filling N to t- 

In Figs. [HJ and [TS] the NRG results of the conduc- 
tances and phase shifts are plotted as functions of ed for 
several fixed values of Se, with St = 0. The values of 
fapcx are chosen such that Se/U = —0.05, —0.1 for Fig. 
[H] and Se/U = 0.05, 0.1 for Fig. [T3] The other parame- 
ters U and r are chosen to be the same as those for Fig. 

El 

We see in (b) and (d) of both Figs. [14] and [15] that the 
series conductance dip emerges in the Kondo plateau at 
three-electron filling N to t — 3.0 even in the presence of 
the small deformations. Particularly, the typical feature 
of the sharp dip clearly remains. The position of the dip, 
however, shifts sensitively to Se, and it moves towards the 
right (left) for Se < (Se > 0). This reflects the behavior 
of the phase shifts at — 0.8(7 < < —0.3(7 in (a) and (c) 



of these figures. Particularly, the odd phase shift S shows 
a clear contrast between the positive and negative Se. 
Correspondingly the kink, which is seen in the even phase 
shift as a step between the values of 2S e /ir = 2.0 and 3.0, 
moves as Se varies. When the gate voltage ed decreases 
from a value at two-electron filling, the third electron 
enters mainly the b\ orbital until ed reaches the SU(4) 
symmetric point, at which point the charge distribution 
changes significantly as seen in Fig. [7] The redistribution 
of the charge is protracted for Se > 0, because the odd 
orbital becomes energetically favorable E^l_ > Eo ■ In 

the opposite case, for < E„ , the redistribution 

takes place earlier than the regular triangle case. As the 
magnitude \Se\ of the deformations increases further, the 
conductance dip moves away from the Kondo plateau and 
disappears. The range of \Se\, in which the dip can be 
seen, will increase with the width of the Kondo plateau. 

Another notable change caused by the deformations 
due to Se is the structure at five electron filling -/V to t — 
5.0. We see in (b) and (d) of both Figs. [HI and [TBI that 
a peak of the series conductance evolves at —1.2(7 < 
e d ^ —1.1(7, as | <5e | increases. Simultaneously, in this 
region, the parallel conductance has a gentle shoulder. 
Furthermore, we can see also in (c) of in these figures 
that a weak step for the five-electron occupancy emerges 
for iVtot in the case of Se = ±0.1(7. 

In the other regions of the filling, namely at iVtot — 4.0 
and iVtot ^ 2.0, the conductances and phase shifts are less 
sensitive to the small perturbations due to Se. Specifi- 
cally, the S = 1 Kondo behavior is robust, although the 
region of ed for the four-electron filling becomes slightly 
narrow. As mentioned in Sec. Ill Bl the Nagaoka state 
is the ground state of the TTQD for AT to t = 4 in the 
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FIG. 16: (Color online) Phase shifts and conductances as 
functions of ed, for (upper panels) t' = l.Qlt and (lower pan- 
els) t' = 0.98i, with U/(2irt) = 1.0, Y/t = 0.12 and e apcx = e d . 



limit of F = 0. We have calculated the cluster eigen- 
states in the presence of the deformations, and found 
that the Nagaoka state remains as the ground state 
near ed/U ~ —1.0 in a wide range of the deformation 
-0.41*7 < Se < 0.87*7 for U/(2nt) = 1.0. This is 
determined mainly by the finite-size energy separation 
AE^ = E { S 4 1 - between the Nagaoka and the 

lowest singlet states, mentioned in Sec. Ill Bl This energy 
scale of the TTQD makes the 5 = 1 Kondo behavior 
rigid also against the deformations. 

We have examined also the effects of small deforma- 
tions caused by the hopping matrix elements St = t' — t, 
choosing the level positions to be uniform Se = 0. In Fig. 
[TBI the conductances and phase shifts are plotted as func- 
tions of td for t' = l.Olf (upper panels), and for t' = 0.98i 
(lower panels) . Besides huge similarities to the above re- 
sults for the case of Se, one notable difference is that no 
pronounced structures evolve at N tot = 5.0 for small St. 
In the cluster limit r = 0, the five-electron charge state 
can become the ground state for finite St, but in a very 
narrow region near ed — —1.16*7, the width of which Aed 
is calculated to be Ae d ^ 0.004*7 for t'/t = 0.98, and 
Ae d ~ 0.002*7 for t'/t = 1.01. This energy scale is much 
smaller than the hybridization strength T = 0.019*7, and 
thus quantitatively the five-electron structure was not 
seen in the present situation. 

In the other regions of the filling, however, the re- 
sults obtained for St > and St < resemble those for 
Se < and Se > 0, respectively. Particularly, the struc- 
tures of the series conductance dip and the phase shifts 
at iVtot — 3.0 are similar qualitatively to those seen in 
the upper panels of Figs. [Ml and [TBI These features are 
determined essentially by the way that the orbital degen- 
eracy between the even and odd orbitals is lifted, namely 
E^l < EP or E^l > E [ o ] . The SU(4) Kondo effect 



is sensitive to relatively small deformations, because it 
relies on the orbital degeneracy at three-electron filling. 
Nevertheless, the sharp structure of the dip does not col- 
lapse immediately under an infinitesimal perturbation, 
and thus a perfect C^ v symmetry is not necessary to ob- 
serve the SU(4) Kondo behavior in the TTQD. 

The S — 1 Kondo behavior at the four-electron fill- 
ing is not sensitive to St, as that in the case of small Se. 
We have confirmed also that the Nagaoka state remains 
as a ground state of the isolated TTQD cluster near 
ed/U ~ —1.0 also in a wide of range of the off-diagonal 
inhomogeneity 0.30 < t'/t < 2.87 for *7/(27rf) = 1.0. The 
5 = 1 Kondo behavior is robust also against the deforma- 
tions by St, and it is mainly due to the finite-size energy 
separation of AE^ . 



VII. SUMMARY 

We have discussed first of all in Sec. Mil that the fixed- 
point Hamiltonian which describes the low-lying energy 
excitation from a local Fermi-liquid ground state has 
an SU(4) symmetry, consisting of the channel and spin 
degrees of freedom, at zero points of the two-terminal 
conductance, or equivalcntly at the points where the 
two phase shifts satisfy the condition 5 e — S = nir for 
n = 0,±1,±2, . . .. This symmetry property holds at low 
energies quite generally for the systems which can be de- 
scribed by the Hamiltonian TL given in Eq. UJ, for any 
numbers of quantized levels Njj, in the case that the sys- 
tem has time-reversal and inversion symmetries. Specif- 
ically, the Hamiltonian 7i does not necessarily have an 
SU(4) symmetry on the global energy scale. 

In the main part of the paper we have discussed the 
Kondo effects taking place in the triangular triple quan- 
tum dot connected to two non- interacting leads, as illus- 
trated in Fig. [TJ We have deduced from the difference 
between the two phase shifts shown in Fig. [5] that the 
two-terminal conductance of this system can have three 
zero points as the level position ed varies. The low-energy 
properties at these zero points can be characterized by 
the SU(4) symmetric local Fermi-liquid theory. Particu- 
larly, the zero point at three-electron filling appears as a 
sharp dip in the middle of the conductance plateau seen 
in Fig. \5\ At this zero point the phase shifts satisfy the 
conditions 5 e — 5 Q = n and S c + S ~ 37r/2, which means 
that S e — 57r/4 and S ~ 7r/4. Furthermore, at this 
filling, the system has an SU(4) symmetry also at high 
energies, which is caused by the four-fold degeneracy in 
the local states of an isolated triangular cluster. These 
symmetry properties at two opposite limits, at low and 
high energies, cause a single-stage screening observed in 
temperature dependence of the entropy shown in Fig. 1111 

For the values of ed away from the dip, the system no 
longer has the low-energy SU(4) symmetry, and the low- 
lying energy excited states have only the SU(2) symmetry 
of the spin. We have confirmed also that the charge dis- 
tribution in the even and odd orbitals, the geometrical 



16 



feature of which is illustrated in Fig. [3l changes signifi- 
cantly at the dip as shown in Fig. [7J The inhomogeneous 
charge distribution in the even-odd basis also causes the 
two-stage Kondo behavior of the S = 1 local moment 
which emerges at four-electron fillings due to Nagaoka 
ferromagnetic mechanism. The partial moment which is 
distributed away from the two leads makes the Kondo 
temperature for the second stage very small. Neverthe- 
less, the Kondo temperature can be raised by increasing 
the coupling T between the TTQD and leads. Therefore, 
in order to observe the second screening stage, T should 
be large. At the same time, however, T has to be smaller 
than AE^ ~ max([//3, t) so that the hybridization ef- 
fects should not smear the typical structures of the 5=1 
Kondo behavior. 

We have studied also the effects of the deformations, 
which break the full triangular symmetry C% v and lift 
the orbital degeneracy in the one-particle excited states. 
Specifically, the deformations which are caused by inho- 
mogeneity in the onsite potentials Se = e ap cx — £d and 
the inter-dot couplings St = t' — t have been examined. 
We found that the series conductance dip does still ex- 
ist stably at three-electron filling, even in the presence 
of infinitesimal deformations. The position of the dip, 
however, shifts sensitively from the middle of the Kondo 
plateau for N to t — 3.0. In the other region of the fillings, 
small deformations do not affect the Kondo behavior sig- 
nificantly, and just cause some gradual and quantitative 
changes. These results support our belief that the var- 
ious Kondo effects in the TTQD could be observed in 
future. 
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APPENDIX A: NAGAOKA STATE IN A 
TRIANGLE 

The eigenvector of the Nagaoka state I^^Lj) for the 
isolated triangle with t' = t and e a pcx = £d can be con- 
structed simply, by filling first the one-particle state of 
k = defined in ([6]) by two electrons, and then filling each 
of the two excited states with k = +2tt/3 and k — —2tt/3 
by a single up-spin electron. This state can be expressed, 
using the even and orbitals described in Eq. (jTSJ) and Fig. 
[3l in the form 

l41i)= a \I) + ai \IT) , (Al) 
|I) = &ULat T oS,|0), |II) = 6t T aLaLaJ,|0>. (A2) 



Here, ao = y/T/3 } a\ = v/2/3, |0) is a vacuum, and 

the eigenvalue is give by £5^1 = —2t+U + Aed- This 
state has a finite total-spin z component of S z = +1, and 
breaks the spin rotational symmetry. In this state, the 
partial moments for even orbitals are given by 

m a>i = (<s> ( sLi\ 2 (4r a *t - 4i a u) 

= \ (1 - H 2 ) • (A3) 

Thus, TO a ,o = 1/3, ™ a ,i = 1/6, and the partial moment 
in the odd orbital b 1 is given by rrib.i = 1/2. Furthermore, 
the occupation number of each orbital takes the value 
(n a ,o) = 4/3, (n a ,i) = 5/3, and (nb,i) = 1. 

In the subspace with Y to t = 4 and S — 1, the matrix 
elements of the interaction Hamiltonian has a diagonal 
form 

{v\H u Aot \v') = UKy , (A4) 

for v, v' = I, II. This can be verified directly using Eq. 
(fT!)]) . For this reason, the wavefunction is determined by 
the noninteracting Hamiltonian, which can be rewritten 
generally in the following form (see also Fig. [3]) 

^dot = e a P cx ?Vo + (e<J - t') n ,i + (ed + t') n M 

-V2t^2 + °>L%*) ■ (A5) 

Therefore, the coefficients ao and ot\ defined in Eq. (|A2|) 
depend crucially on the potential difference between the 
a and a 1 orbitals, which correspond to the first and sec- 
ond terms in the right-hand side of Eq. (|A5[) . This po- 
tential profile in the even-odd basis play an important 
role on the two-stage Kondo effects, which take place in 
the case that the two leads are connected to the a l and 
b 1 orbitals in the way that is shown in Fig. [3l 

APPENDIX B: MANY-BODY PHASE SHIFTS 

The phase shifts for interacting electrons can be de- 
fined, using the Green's function 

Gij(i0J n ) = - j^dr (T T d ic (r) dJ^O)) , (Bl) 

where (3 = 1/T, dj a (r) = e TH d 3<y er TH , and (O) = 
Tr \e~ m O\ /Tre-P n , The retarded and advanced 
Green's functions correspond to the analytic continua- 
tions Gfj(uj) = Gij(uj ± i + ). The interaction H% ot be- 
tween electrons in the dots at 1 < i < (Nn = 3 for 
the triple dot) causes the self-energy correction Hij(z), 
and the Dyson equation is given by 

N D N D 

Gijiz) = G%{z) + £ E G °"'^ ■ 

i'=l j' = l 

(B2) 
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Here, G^j(z) is the non-interacting Green's function cor- 
responding to H = ^d ot + U mh[ + U\ cad . 

At zero temperature T = 0, the series g s and parallel 
<7 P conductances can be expressed in terms of the Green's 
function at the Fermi level to = 0, in the case the system 
has a time-reversal symmetr y 19 ' 23 



9 P 



^-4r fl r L |G+ Dl (0)| 2 

h 



g+(o)| |g+ d7Vd (o)| 



2e2 
h 



2T L T R |G+ Dl (0)| 



-r L imG+(o)-r fl im g+ dJVd (o) 



(B3) 

(B4) 
(B5) 



Note that the contributions of the vertex correction do 
not appear here due to the property that the imaginary 
part of the self-energy vanishes ImE^(O) = at T = 
and oj = 0J& Therefore the value of G^(0) appearing in 
the above equations can be reproduced correctly by a free 
quasi-particle Hamiltonian Ti. s = Hdot + 7~tmix + ^load, 



H 



dot 



iVi, 
i,j=l a 



(B6) 



(|B3 |) -([B8 ]) . One can state immediately from Eq. (|B8[) 
that the series conductance g s becomes zero in the case 
k-c = K , or equally for S c — 6 = nit for n = 0, ±1, ±2, . . .. 



APPENDIX C: NRG APPROACH 

In the NRG approach the non-interacting leads are 
transformed into the chains carrying out the logarith- 
mic discretization with the parameter A. and a se- 
quence of the Hamiltonian Hn in the following form is 
introducedr^iS 3 - 

H N = A (JV ~ 1)/2 [H° ot + H^ot + H mix + H^j , (CI) 



a 



Ra J ' 



(C2) 



{N) _ 1 + 1/A 

^lcad - U o 



N-l 

£ EE^" B/2 



v=L,R a n=0 



X ( fn+l,va- fn,va + fn,vcr fn+l,va ) ■ (^3) 



Here, — iy = — ty + e^i tfy+Re (0) can be regarded as 
the hopping matrix element for free quasi-particlesi 19 ' 55 
Furthermore, in the case the system has also the in- 
version symmetry = r# (= T), the Green's functions 
take the form 



Gti(0) = G% DND (0) =-L 



1 



k c + i 
1 

k + i 



K + I 
1 



(B7) 
(B8) 



Here, k c = — cot S c and k q — — cot (5 D include all the 
many-body corrections through the real part of the self- 
energy Re £+(0)^ Equations (|T0 ]) -(P j) follow from Eqs. 



Here, D is the half-width of the conduction band, and 
the other parameters are given b y 52 ' 53 



2DTAt 



AA = i^logA, (C4) 



2 1-1/A 



1 - 1/A 



n.+ l 



y/1 - 1/A 2 "+Vl " 1/A 



2n+3 



(C5) 



The factor A\ represents a correction to the discrete 
model for comparison with results in the continuum limit 
A — > l . 52 ' 56 We have carried out the iterative diagonal- 
ization of Hn using the even-odd basis defined in Sec. 

mm 
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